Transition to Long Range Magnetic Order in the Highly Frustrated 
Insulating Pyrochlore Antiferromagnet Gd2Ti207 
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Experimental evidence from measurements of the a.c. and d.c. susceptibility, and heat capacity 
data show that the pyrochlore structure oxide, Gd2Ti207, exhibits short range order that starts 
developing at 30K, as well as long range magnetic order at T ~ IK. The Curie- Weiss temperature, 
^cw = -9.6K, is largely due to exchange interactions. Deviations from the Curie- Weiss law occur 
below ~10K while magnetic heat capacity contributions are found at temperatures above 20K. A 
sharp maximum in the heat capacity at Tc — 0.97K signals a transition to a long range ordered state, 
with the magnetic specific accounting for only ~50% of the magnetic entropy. The heat capacity 
above the phase transition can be modeled by assuming that a distribution of random fields acts 
on the **S7/2 ground state for Gd^^. There is no frequency dependence to the a.c. susceptibility 
in either the short range or long range ordered regimes, hence suggesting the absence of any spin- 
glassy behavior. Mean field theoretical calculations show that no long range ordered ground state 
exists for the conditions of nearest-neighbor antiferromagnetic exchange and long range dipolar 
couplings. At the mean-field level, long range order at various commensurate or incommensurate 
wave vectors is found only upon inclusion of exchange interactions beyond nearest-neighbor exchange 
and dipolar coupling. The properties of Gd2Ti2 07 are compared with other geometrically frustrated 
antiferromagnets such as the Gd3Ga50i2 gadolinium gallium garnet, RE2Ti207 pyrochlores where 
RE = Tb, Ho and Tm, and Heisenberg-type pyrochlore such as Y2M02O7, Tb2Mo207, and spinels 
such as ZnFe204 
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I. INTRODUCTION 

There has been in the past eight years an enormous 
amount of theoretical and experimental activity devoted 
to the study of highly geometrically frustrated antiferro- 
magnetic materials M . The main reason for this interest 
stems from the suggestion that a high degree of frustra- 
tion can induce sufficiently large zero-temperature quan- 
tum spin fluctuations as to destroy long range Neel order 
even in three dimensions. This could give rise to new 
exotic and intrinsically quantum mechanical magnetic 
ground states such as dimerized ground states, "spin ne- 
matics", or fully disordered states (e.g. RVB-like) with 
no broken spin or lattice symmeties 1^-^!. Frustration 
arises when a magnetic system cannot minimize its total 
classical ground-state energy by minimizing the bond en- 
ergy of each spin-spin interaction individually B. This 
is the case, for example, in systems where antiferromag- 
netically coupled spins reside on a network made of basic 
units such as triangles or tetrahedra. On a triangular pla- 
quette, vector (i.e. XY or Heisenberg) spins can manage 
the frustration better than Ising moments by adopting 
a noncoUinear structure with the spins making an an- 
gle of 120° from each other. Triangular or tetrahedral 
units can be put together to form a regular lattice such 
that they are either edge sharing or corner sharing. For 



example, the space-filling arrangements of edge-sharing 
triangles and tetrahedra form the well-known triangu- 
lar and face-centered cubic lattices in two and three di- 
mensions, respectively. In two dimensions, a network of 
corner-sharing triangles forms the kagome lattice P^pl. 
In three dimensions, a lattice of corner-sharing tetrahe- 
dra forms the structures found in spinels, Lave phases, 
and pyrochlore crystals [0,pP], while corner-sharing tri- 
angles give the familiar garnets p-p^ . 

Among highly frustrated antiferromagnets, the three- 
dimensional pyrochlore lattice of corner-sharing tetrahe- 
dra is aparticularly interesting system (see Fig. 1) . The- 
ory 0,|l^,Q and Monte Carlo simulations ||lj,|l^ show 
that classical Heisenberg spins residing on the vertices of 
the pyrochlore lattice and interacting only via nearest- 
neighbor antiferromagnetic exchange do not show a tran- 
sition to long range magnetic order at nonzero temper- 
ature. This is very different from the, also frustrated, 
classical nearest-neighbor fee Heisenberg antiferromagnet 
where long range order occurs at finite-temperature via 
a 1st order transition driven by thermally-induced order- 
by-disorder |lq-p!8|| . Because of their failure to order even 
at the classical level, the high-frustration present in py- 
rochlore antiferromagnets would appear to make these 
systems excellent candidates to search for novel three- 
dimensional quantum disordered magnetic ground states. 



Indeed, recent numerical calculations suggest that the 
quantum S = 1/2 pyrochlore Heisenberg antiferromag- 
net may be a quantum spin-liquid Q]. 

In all real systems there exist perturbations, {H'}, be- 
yond the nearest-neighbor Heisenberg Hamiltonian such 
as further than nearest-neighbor exchange, single-ion and 
exchange anisotropy, and magnetic dipolar couplings. 
For a classical system, one generally expects that such 
perturbations will select a unique classical ground state 
to which a transition at nonzero temperature can oc- 
cur. It is also possible that the energetic perturbations 
{H'} can sufficiently reduce the classical degeneracy such 
that "additional" order-by-disorder via thermal ||7|, p^ , p0[ 
and/or weak quantum fluctuations ||21|], such as occur 
in the fee antiferromagnet [|l6|-|l8|, can "complete" the 
ground-state selection and give rise to a transition to con- 
ventional long range Neel order at nonzero temperature. 
However, in presence of quantum fluctuations (i.e. fi- 
nite S spin value), one expects that for small spin value 
S and/or sufficientl y w eak {H'}, a quantum disordered 
phase may occur |Q,^ . 

What is perhaps one of the most interesting issues in 
geometrically frustrated antiferromagnet systems, is that 
a large number of highly frustrated antiferromagnetic 
insulators exhibit spin-glass behavior rather spin-liquid 
behavior. For example, in a number of antiferromag- 
netic pyrochlore oxydes, such as Y2M02O7 |23-p4] and 
Tb2Mo207 |2^j2^, spin-glass behavior is observed simi- 
lar to what is seen in conventional randomly disordered 
and frustrated spin glass materials |26|| , even though the 
measured disorder level is immeasurably small in the py- 
rochlores p7| , p8| . 

Pyrochlores oxides, of which two examples have just 
been cited, present a number of opportunities for study- 
ing geometrically frustrated antiferromagnets. In this 
structure both cation sites in Fd3m, the 16c site nor- 
mally occuped by a transition element and 16d normally 
occupied by a lanthanide, have the same "pyrochlore" 
topology, i.e., the three dimensional network of corner 
sharing tetrahedra shown in Figure 1. Thus materials 
exist with only the 16c site magnetic (Y2M02O7) p2|-p4[, 
with both sites magnetic (Tb2Mo207) ||2|,|2|, and with 
only the 16d site magnetic, the RE2Ti2 07 materials for 
example. This latter series has attracted attention very 
recently with studies for RE= Tm, Ho and Tb [||-||. 
The RE2Ti207 materials are generally quite well-ordered 
crystallographically, with oxygen non-stoichiometry and 
16c/16d site cation admixing at or below the limit of 
detection by neutron diffraction p8[ . 

Gd2Ti207 has not yet been studied in detail and there 
exist some compelling reasons to do so [Q . Previous re- 
ports for Gd2Ti2 07 indicate no long range order down 
to IK [|||. The Gd3+ ion (if) is spin-only with a 
^5*7/2 ground state and thus, crystal field splittings and 
anisotropy, which play a large role in the properties of 
the aforementioned Tm'^^, Ho'^"'" and Tb'^"'" materials, 
as will be discussed later, are expected to be relatively 
unimportant. Gd2Ti2 07 then, should be an excellent 



approximation to a classical Heisenberg antiferromag- 
netic system with dipole-dipole interactions as leading 
perturbations {H'}. In addition it is important to com- 
pare this material to the Gd3Ga50i2, gadolinium gallium 
garnet (GGG), where the Gd'^^ ions reside on a three- 
dimensional sublattice of corner-sharing triangles p-fi2[ . 
GGG has been found to possess a very unusual set of 
thermodynamic properties with anomalous specific heat 
behavior, spin glass magnetic properties and no true long 
range order |y,^Q], but incommensurate short range or- 
der developing at very low temperatures [0| . Interesting 
properties are also observed for applied magnetic fields 
in the range [0.1 — 0.7] Tesla [H. same time Monte 
Carlo simulations for GGG have found some intrinsic 
(e.g. disorder-free) glassy behavior with no evidence for 
the development of short range order 10,0]. Hence, one 
further motivation for studying Gd2Ti2 07 was to hope 
bridge a gap between the peculiar behavior of GGG and 
the transition metal pyrochlores, such as Y2M02O7 and 
Tb2Mo207 which show spin-glass behavior and, hence, 
possibly indirectly gain some insight on the thermody- 
namic behavior of both GGG and the insulating antifer- 
romagnetic Heisenberg pyrochlore systems. 

In this work a detailed study of Gd2Ti2 07 has been 
carried out including both a.c. and d.c. susceptibility 
and heat capacity studies. To complement the experi- 
mental work, results from mean-field theoretical calcu- 
lations are presented which take into account exchange 
and dipolar interactions. 



II. EXPERIMENTAL METHOD 

A. Sample Preparation 

A polycrystalline sample of Gd2Ti2 07 was prepared by 
high temperature solid state reaction. Starting materials, 
Gd203 and Ti02, were taken in stoichiometric propor- 
tions and mixed thoroughly. The mixture was pressed 
into pellets and heated in an alumina crucible at 1400C 
in air for 12 hours. The powder x-ray diffraction pattern 
of the sample obtained using a Guinier-Hagg camera in- 
dicate that the sample formed is single phase with the 
cubic pyrochlore structure. The size of the conventional 
cubic unit cell is gq = 10. 184(1) A. 



B. DC & AC Magnetic Susceptibility Measurements 

The DC magnetic susceptibility, x, was measured using 
a SQUID magnetometer (Quantum Design, San Diego) 
in the temperature range 2— 300K. The AC susceptibil- 
ity, xaCi was measured at different frequencies by the 
mutual inductance method. The primary coil of the mu- 
tual inductor is energized by a frequency generator (DS 
335, Stanford Research Systems) and the output across 



the two identical secondary coils, wound in opposite di- 
rections, was measured using a lock-in-amplifier (SR-830 
DSP, Stanford Research Systems). The sample suscepti- 
bility was determined from the difference in the outputs 
with the sample in the middle of the top secondary coil 
and without the sample. The cryostat used for the tem- 
perature variation is described in the section below. 



C. Specific Heat Measurements 

The specific heat of the sample in the form of a pel- 
let («100 mg) was measured in the temperature range 
0.6— 35K using a quasiadiabatic calorimeter and a com- 
mercial Heliox sorption pumped '^He cryostat supplied 
by Oxford Instruments. The sample was mounted on a 
thin sapphire plate with apiezon for better thermal con- 
tact. Underneath the sapphire plate a strain gauge heater 
and a RUO2 temperature sensor were attached with G- 
E varnish. The temperature of the calorimetric cell was 
controlled from the ^He pot on the Heliox. The sample 
temperature was measured using an LR-700 AC resis- 
tance bridge at a frequency of 16Hz. The specific heat of 
the sample was obtained by subtracting the contribution 
of the addendum, measured separately, from the total 
measured heat capacity. 



III. EXPERIMENTAL RESULTS 

The DC susceptibility. Fig. 2a, x, measured at an 
applied field of 0.01 Tesla vs temperature is found to 
obey the Curie- Weiss behavior in the range lOK— 300K. 
An effective magnetic moment of 7.7/iB/Gd'^^ obtained 
from the Curie- Weiss fit is close to the expected value 
of 7.94^B/Gd"^^ for the free ion, ^S7/2, and a param- 
agnetic Curie temperature, 6*0-^1 of -9.6(3)K indicates 
antiferromagnetic interactions between the Gd'^"'' spins. 
It is worth noting that x starts deviating at a temper- 
ature of the order of 0cw as it ought to be for a "con- 
ventional" system undergoing a transition to long range 
order. That ^cw is predominantly due to exchange in- 
teractions as opposed to crystal field effects is confirmed 
by measurements on the magnetically diluted system 
(Gdo.o2Yo.98)2Ti207, for which 9cw is much reduced and 
of the order of ^ — 0.9K (Fig. 2b). The absence of any 
magnetic ordering down to 2K in the concentrated sys- 
tem, even though 6'cw is about five times larger than this 
temperature, suggests the presence of important mag- 
netic frustration inhibiting the occurrence of magnetic 
long range order. 

In search of a possible magnetic ordering below 2K, 
AC susceptibility, xaCj was measured down to 0.3K. The 
temperature variation of xac for different frequencies. 
Fig. 3, exhibits two features, a broad peak centered at 
about 2K and a sharp down turn below about IK, the 



latter possibly signaling a transition to long range antifer- 
romagnetic order. XAc('^) appears to be independent of 
frequency which would seem to rule out a spin glass state, 
as opposed to what has been found in other pyrochlore 
oxydes such as Y2M02O7 ^^, Tb2Mo207 ^,|| and 
the frustrated Gd3Ga50i2 garnet §,|o|jl|. 

The specific heat, Cp, as a function of temperature is 
shown in Fig. 4. There is a broad peak centered around 
2K and a very sharp peak slightly below IK indicating 
the presence of short range correlations and, in agree- 
ment with the AC susceptibility data, the development 
of long range magnetic order via a sharp transition at IK. 
The solid line corresponds to the estimated lattice specific 
heat, Ci, of Gd2Ti207 determined by scahng the specific 
heat for Y2Ti207, which is insulating, non-magnetic, and 
isostructural to Gd2Ti207. The magnetic specific heat. 
Cm 1 was obtained by subtracting C'l from Cp and its tem- 
perature variation is shown in Fig. 5. 

The Gd'^+ ion has an isotropic spin of S* = 7/2 with no 
orbital magnetic moment contribution and the degener- 
acy of the 8— level ground state cannot be lifted by the 
crystal electric field beyond a fraction of a Kelvin. The 
presence of Cm up to about 30K clearly indicates that 
the ground state degeneracy is lifted by magnetic inter- 
actions. 

The magnetic entropy, Sm, was obtained by extrapo- 
lating the Cm/T behavior to OK and numerically inte- 
grating it vs temperature. The total magnetic entropy is 
33.8 J mol-lK-1 which is close to the expected 2Rln(8) 
=34.6 Jmol^^K"^ for an 5" = 7/2 system. The entropy 
recovered at the long range order temperature is about 
50% of the total value which indicates that a sizeable 
fraction of the entropy is due to the short range correla- 
tions present above T = IK. 

An attempt was made to fit the Cm data above IK. 
A zeroth order model consisting of a simple Schottky 
anomaly based on a splitting scheme of 8 equally spaced 
discrete levels, i.e., assuming a unique value for the in- 
ternal magnetic field at each Gd'^+ site, gives a poor fit, 
as might be expected. A much better fit is obtained by 
assuming a continuous range of energy level splittings 
with a truncated Gaussian distribution. The probability 
distribution is normalized such that the area under the 
curve is unity. The resulting fit is shown in Fig. 5. This 
model is equivalent to assuming a distribution of inter- 
nal magnetic fields, i.e., an array of random fields as ap- 
propriate to a thermal regime dominated by quasi-static 
short range magnetic order. A similar approach has been 
used before to model the specific heat anomaly due to 
the Gd sublattice in Gd2Mo207 [g^l- I^ this pyrochlore 
structure material the Mo^+ sublattice undergoes a spin 
glass type of order at 60K while the Gd'^^ specific heat 
contribution is also a broad peak centered about 9K. 

In summary, experimental results obtained from AC 
and DC susceptibility measurements as well as specific 
heat measurements reported in Section III give strong 
compelling evidence for a single sharp transition to a 
long range ordered state at Tc = 0.97K preceeded by 



a short range ordered regime which extends to approxi- 
mately 30Tc ~ 36'CW. The heat capacity in this regime 
can be modelled in terms of a distribution of random ex- 
change fields acting on the ^87/2 ground state of Gd'^+. 



IV. MEAN-FIELD THEORY 
A. Model and Method 

Our aim in this section is to determine, within mean- 
field theory, what the expected magnetic properties and 
type of magnetic ordered phase(s) for a classical spin 
model of Gd2Ti207 are. 

We first consider the following classical spin Hamilto- 
nian for Gd2Ti207: 



pyrochlore antiferromagnets ||I^. We consider the 
mean-field order parameters, B(ri) at site r^, in terms 
of Fourier components. The pyrochlore lattice is a 
non-Bravais lattice, and we use a rhombohedral ba- 
sis where there are four atoms per unit cell located 
at (0, 0, 0), (1/4, 1/4, 0), (1/4, 0, 1/4), and (0, 1/4, 1/4) in 
units of the conventional cubic unit cell. We relabel the 
spins, S(ri), in terms of unit cell coordinates, and a sub- 
lattice index within the unit cell, and take advantage of 
the translational symmetry of the lattice, and expand the 
order parameters B(ri) in terms of Fourier components. 
In this case B"(ri) on the a'th sublattice site of the unit 
cell located at r, can be written as 



B'^(rO=^B'^(q)exp(zq.O 
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The first term is the isotropic Heisenberg exchange in- 
teraction, and the second term is the dipolar coupling 
between the Gd magnetic moments. For the open py- 
rochlore lattice structure, we expect very small sec- 
ond and further nearest-neighbor exchange coupling, 
Jn>2, compared to the nearest-neighbor Ji (J„>2 < 
O.OSJi) [Q. Hence, we first consider the case where the 
sum in the first (exchange) term of Eq.(4.1) above is re- 
stricted to the nearest-neighbor exchange Ji only |3^. 

Gd'^^ has a spin S = 7/2, which gives an effective 
dipole moment of /x(Gd^+) = g^B\/S{S + 1) = 7.94^b, 
with 5 = 2, in good agreement with the Curie constant 
determined in Section III. This gives an estimate for the 
nearest- neighbor strength of the dipole-dipole interaction 
Ddd — 63/iB/io/(47rrfj„), where /xq is the magnetic per- 
meability. With a Gd^+ at r = (0,0,0) and a nearest- 
neighbor at r„„ = (a/4, a/4, 0), where a = 10.184A 
is the size of the conventional cubic unit cell, we find 
Ddd ~ 0.84K. 

An estimate of the nearest-neighbor exchange Ji can 
be found from the measured Curie- Weiss temperature 
(see below). We have ^cw ^ — 9.6K. This gives for 
the effective classical nearest-neighbor exchange, Jf = 
JiS{S -t- 1) - -4.8K using 6'cw = zJiS{S + l)/3, where 
z = 6 is the number of nearest- neighbor. We hence- 
forth use a classical approximation of Eq. (2.1) above, 
where we use unit length vectors Si, and replace ^i 
by SiiDddY , and express Vij in units of the nearest- 
neighbor distance. We set Ji — — 4.8K and a strength of 
0.84K for the nearest-neighbor dipole coupling. Below, 
we take Ddd/Ji = 0.2 ^^. Hence, unhke in the transi- 
tion metal pyrochlores, dipole-dipole interactions is the 
major perturbation at play beyond the nearest-neighbor 
Heisenberg exchange coupling in Gd2Ti2 07. 

We now proceed along the lines of Reimers, Berlin- 
skly and Shi in their mean-field study of Heisenberg 



The spin-spin interaction matrix, Jai3i\'<^ij\), including 
both exchange and dipolar interactions, reads: 



JS^fli\^^j\) = JlSal3Sr,,,r„„ + Ddd ' 
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where 6af3 is the Kronecker delta, and a and /? refer to 
the x,y,z cartesian components of Si and r^j". rf^"^^ de- 
notes the a components of the interspin vector r^ that 
connects spin S" to spin S^. We write j7a/3 in terms of 
its Fourier components as 

q 

where N is the number of unit cells with 4 spins per unit 
cell. 

The quadratic part of the mean-field free-energy, F^'^' , 
then becomes ||k| : 

F(2)(T)/7V = i Yl ^"(q) {^TSabSo.^ - J:^(q)} i?^(-q) 
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Diagonalizing F^'^' (T) requires transforming to normal 
modes of the system 
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where $fl(q) are the eigenmodes, and U{q) is the unitary 
matrix that diagonalizes J{c\) in the spin— sublattice 
space, with eigenvalues A(q) 
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Henceforth we will use the convention that indices (ab) 
label sublattices, that indices {ijk) label the normal 
modes, and that (a/37) label spin components. We ex- 
press F^^) (T) in terms of normal modes 

^^'VA^ = jEEE'^;(q)*7(-q){3r-^;(q)} 



q I 7 
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The first ordered state of the system occurs at the tem- 
perature 



T, = imax{AJ.(q)} 

O q.4,a 



(4.9) 



where maxq ^ Q{A^(q)} indicates a global maximum of 
the spectrum of A„ (q) for all q. 

Let us briefly explain how we procede using the above 
set of equations to determine the "soft mode(s)" of the 
system at Tc- Firstly, the Fourier transform of J^'sil^ijl) 
is calculated using Eq. (4.3) ^]. For the rhombohedral 
basis used above, the space is of dimension Vs ® T>s\, 
where the spin-component subspace, Vs, is of dimension 
3x3 and the sublattice subspace, 2?sb is of dimension 
4x4. The eigenvalues, {AJj(q)}, and eigenvectors <&a(q) 
are determined by reshaping J7^«(q) into a 12 x 12 array. 
The pyrochlore lattice has a symmetry of inversion with 
respect to a lattice point and this imphes that J^^(q) 
is real and symmetric. The eigenvalues and eigenvectors 
are found using a standard numerical packages for eigen 
problems of real symmetric matrices. 



B. Results 

For Ddd = 0, we recover the results of Ref. [^. Before 
we present the results with the dipolar interactions, we 
review what the mean-field results found for the isotropic 
pyrochlore class Heisenberg antiferromagnet depending 
on the values of the second, J2 and third, J3 nearest- 
neighbor exchange couplings are fl^. For J2 = J3 = 
there are two dispersionless unstable or critical modes 
throughout the Brillouin zone. There are therefore no se- 
lected wavevector for long range order. Numerical work 
has shown that no long range order occurs at nonzero 
temperature in the nearest-neighbor classical Heisenberg 
pyrochlore antiferromagnet |14 1^. For J3 = 0, ferro- 
magnetic J2 > gives rise to an ordering at an in- 
commensurate wavevector, while for antiferromagnetic 
J2 < 0, the system orders at q* = 0. For J2 — 0, 
and ferromagnetic J3 > 0, the system also orders at 
q* = 0, while there are dispersionless (degeneracy lines) 
along certain symmetry directions for J2 — and anti- 
ferromagnetic J3 < 0. In the overall parameter space 
{J2/Ji,J3/Ji}, fong range order is always expected to 



occur at nonzero temperature within mean-field theory, 
except for J2 = and antiferromagnetic J3 < (Fig. 6 
in Ref. |l|]). 

We now consider the case where Ddd/ Ji — 0.2, and 
first set J2 = J3 = 0. Naively, one might have thought 
that (i) the long range and (ii) anisotropic nature of the 
dipolar interactions would lift all macroscopic ground 
state degeneracies that occur in the isotropic nearest- 
neighbor (Ji < 0) Heisenberg antiferromagnet and give 
rise to a unique selected wavevector q*, at which long 
range order would occur Q. This is not the case. We 
find that the largest eigenvalue AJj(q) that controls the 
mean-field Tc (Eq. 4.9) is dispersionless along the star of 
the [111] direction in the cubic basis (Fig. 6). The figure 
shows Amax(q) as a function of qi in the [110] and q2 
in the [001], where Aniax(q) is the largest eigenvalue of 
A„(q) at a given q. Hence, no-long range order is to be 
expected in this system within the mean-field approxi- 
mation [}4l|| . In this context, in it interesting to note that 
the combined long range dipolar and RKKY interactions 
in the problem of nuclear magnetism in Cu and Ag do 
not lead either to a full selection of a unique classical long 
range ordered state below the mean-field Tc ^^. Such 
"degeneration lines" as found in the present system also 
occur in other frustrated systems such as the nearest- 
neighbor Heisenberg fee antiferromagnet where there are 
degeneration lines along the 7r/a(l, q, 0) direction [Qq nSJ. 
Degeneration (spiral) lines also occur in the more com- 
plicated case of the rhombohedral antiferromagne t [ 43[ . 
By analogy with the work on the frustrated fee [^ [l8l42[ 
and rhombohedral [E3[ antiferromagnets, we expect that 
for deneration lines (as opposed to degeneration zone as 
in the case of the nearest-neighbor pyrochlore antiferro- 
magnet [ p^ , p^ ). Thermal and/or quantum fluctuations 
will restore long range order at finite temperature via a 
process of order-by-disorder. Work in that direction is in 
progress and will be reported elsewhere ||39|| . 

We find that for either nonzero ferromagnetic or an- 
tiferromagnetic, IJ2I -C Ddd and/or IJ3I -C Ddd, that 
the line-degeneracy along the [111] direction is lifted and 
that a specific value q* along that direction is picked- 
up, giving rise to an absolute maximum of A^(q). For 
sufficiently large | J2I and/or | J3I compared to Ddd, a dif- 
ferent selected wavevector direction is chosen as found by 
Reimers et al. [Q, except that here, there is no degen- 
eration line occuring for J2 = and J3 < 0, as found in 
Ref. H when Ddd ¥^ 0, Ddd < | J2I and Ddd<^\h\. In 
other words, all non-global degeneracies are lifted in the 
case where Ddd 7^ 0, J2 7^ 0, and J3 7^ 0. 

In summary, we would expect that long range or- 
der should occur in Gd2Ti207, either via an order-by- 
disorder mechanism, or via energetic selection of an or- 
dering wavevector via superexchange couplings beyond 
J I and dipolar interactions, T)dd- 



V. DISCUSSION 

It is useful to compare Gd2Ti207 with related systems 
such as the remaining RE2Ti207 materials, Gd3Ga50i2 
(GGG), Gd203, cubic-Gd203 (C-GdzOs) and also 
GdAlOa for reasons which should soon become clear. 

C-GdzOs crystallizes in the so-called bixbyite struc- 
ture, Ia3, with two distinct crystallographic sites. How- 
ever, the sublattice of the two sites taken together is an 
excellent approximation to an fee lattice. Its properties 
are as follows. Curie- Weiss behavior is observed with 
^cw = -17K but magnetic order of undetermined range 
is not found down to 1.6K |44). In fact the description of 
the neutron diffraction data for C-GdzOs bears a strik- 
ing resemblance to that for GGG [Q . The heat capacity 
from 1.4K to 18K shows only a broad "Schottky" type 
anomaly peaked near 2K with magnetic contributions ex- 
tending to about 20K B3]. Neutron diffraction data in 
the form of diffuse scattering confirm that the short range 
magnetic correlations do extend to at least 20K [Q . It is 
clear that C-GdzOa should be reconsidered as a geomet- 
rically frustrated antiferromagnet material although the 
original interpretation of its properties was not presented 
in those terms. In particular it is important to determine 
whether C-GdzOa does indeed undergo true long range 
order as seems to be the case for GdzTizOy and if the 
susceptibility is frequency dependent. 

GGG has been recognized as a geometrically frustrated 
antiferromagnet since 1979 0. Its properties are an 
amalgam of those of GdzTizOy and C-GdzOs scaled to 
lower energies. For example 6'cw = -2.3K and the maxi- 
mum in the heat capacity occurs at 0.8K. True long range 
order is not established down to 43mK but incommensu- 
rate short range order on a ~ 100 A length scale is found 
at the lowest temperature studied |1^]. The extended 
short range order starts developing rapidly at '^ 140 mK, 
which is close to the temperature at which the nonlinear 
coefficient Xsi^) peaks as found in Ref. [Q. Coexis- 
tence of spin glass and long range order in the form of a 
frequency dependent susceptibility is found, which con- 
trasts sharply with the lack thereof in GdzTizOy. At this 
point it is useful to point out that GdzTizOr, GGG and 
C-GdzOa represent anomalies in the context of magnetic 
Gd"^+ oxides. For example monoclinic GdzOs in which 
the Gd sublattice is not fee, unlike C-GdzOs, is a normal 
antiferromagnet with Tc— 3.8K and greater than 90% 
entropy removal below Tc p9]. In GdAlOa the Gd ions 
are on a simple cubic lattice and 6'cw — -4.8K, Tc = 
3.8K and nearly 100% of the entropy is removed below 

At present the strong contrast in behavior between 
GdzTizOT and GGG remains unexplained. On the sub- 
ject of why GdzTizOy orders and GGG does not, one can 
only speculate. For example, for the very specific topol- 
ogy of the garnet lattice, there may remain degenerate 
or quasi degenerate dispersion lines or surfaces of zero 
mode in g-space which survive even upon the inclusion 



of perturbations such as dipolar and or higher than first 
neighbor exchange interactions. Thus the selection of an 
ordering wave vector, q*, for GGG may be much less 
robust than for the much different topology presented 
by GdzTizOy for perturbations {H'} of similar order of 
magnitude P9|] . 

Returning to the REzTizOy series, as mentioned, 
GdzTizOy offers the opportunity to study a system in 
which the crystal field and anisotropy perturbations are 
minimized. As this is certainly not the case for RE= Tb, 
Ho and Tm, some comment on the symmetry of the local 
environment at the RE site is in order. The 16d rare- 
earth site is coordinated by two sets of oxygen atoms, 
six 01(48f) and two 02 (8b), giving eight-fold coordina- 
tion overall. It is important to note that the RE site 
symmetry is strongly distorted from cubic, which would 
imply eight equal RE-0 distances (for RE=Gd the sum 
of the ionic radii give 2.42 A) and 0-RE-O angles of 
70.5°, 109.5° and 180°. In REzTizOr, the six 01 atoms 
form a puckered ring about the RE (Gd-0 distance of 
2.55A) and the two 02 atoms a hnear 02-RE-02 unit 
oriented normal to the mean plane of the puckered ring 
with extremely short RE-02 distances(Gd-0 distance is 
2.21A). This Gd-0 distance is among the shortest, if not 
the shortest, such distance known in Gd oxide chemistry 
and implies a very strong interaction. This observation 
suggests that a crystal field of axial symmetry might be 
an even better approximation than cubic. The 02-Gd- 
02 angle is of course 180° while the 01-Gd-Ol angles 
are 62.7°, 117° and 180°. and the Ol-Gd-02 angles are 
81° and 100°. Thus, it is best,when thinking about the 
RE site crystal field, to consider the true symmetry, - 
3m(D3d), rather than relying on cubic or axial approxi- 
mations. 

The known situation with respect to the presence or 
absence of long range order in the REzTizOy pyrochlores 
is summarized in Table 1. 

One obvious correlation is that those REzTizOy py- 
rochlores which contain a Kramers (odd electron) ion 
nearly always show long range order (Dy being the excep- 
tion) while those with a non-Kramers (even electron) ion 
do not. A "zeroth order" interpretation of the trends 
in Table 1. is then, that the action of the relatively 
low symmetry crystal field induces a true singlet ground 
state in the non-Kramers ions and this is the explana- 
tion of the absence of long range order (LRO). There is 
good evidence that such is the case for RE=Tm from 
a combination of susceptibility |29y30f| , inelastic neutron 
scattering [gO|, and crystal field calculations (using the 
correct -3m symmetry) ^^l- Experimentally |3y], the 
singlet state is well-separated by 120K from the nearest 
excited state which is in remarkable agreement with the 
aforementioned crystal field calculations which predict 
118K |§. 

The other two non-Kramers ions are not so simple. For 
HozTizOy the ground state is thought to be an Ising dou- 
blet p3| , in agreement with crystal-field calculations Ew , 
and the nearest neighbor exchange is weakly ferromag- 



netic. Here it has been argued that the strong Ising- 
Uke single ion anisotropy along the [111] frustrates the 
development of long range ferromagnetic order pi| , p2[ . 
However, recent studies suggest a more complex picture 
where dipolar interactions competing with antiferromag- 
netic exchange is responsible for the behavior observed in 
Ho2Ti207 pM. This material also exhibits spin dynamics 
and spin freezing reminiscent of the disorder-free, intrinsic 
glassy behavior exhibited by the "ice model" [ p]j32| 
with an exponential decrease of the spin lattice relax- 
ation rate suggestive of Orbach processes [p3. In con- 
trast Gd2Ti207 exhibits no apparent dynamics or spin- 
glassiness at any temperature even above Tc- 

A detailed study of Tb9Ti207 will be described in a 
subsequent publication [^ . The salient facts are that the 
Tb'^^ ground state also appears to be a doublet but not so 
well isolated from several other levels within 15K— lOOK. 
The exchange interactions are relatively strongly antifer- 
romagnetic, comparable to Gd2Ti207, and short range 
magnetic correlations persist up to at least 30K, also sim- 
ilar to what is found in Gd2Ti207. Tb2Ti207 does not 
order down to 70 mK pi] . The lack of LRO in this sys- 
tem is difficult to understand pq ]. Indeed, as argued in 
refs [ pO| , ^ a nearest-neighbor Heisenberg antiferromag- 
net with a [111] easy axis is a trivial problem with an 
effectively non-frustrated and unique (two-fold Ising-like 
globally degenerate) classical ground state, and should 
therefore show a phase transition at nonzero tempera- 
ture in the limit of sufficiently strong crystal-field level 
splitting compared to the superexchange J. 

From the above discussion we can conclude that each 
RE2Ti207 material presents its own special set of cir- 
cumstances where details of the finely tuned relative 
strength of crystal field parameters, exchange and dipo- 
lar couplings play a crucial role, and a blanket expla- 
nation for the apparent systematics of Table 1 will not 
be found. It is worth noting an interesting paradox. 
Gd2Ti207 represents the case for which some of the per- 
turbations which might be thought to aid in the selec- 
tion of an unique ground state, i.e., crystal fields and 
anisotropy, are largely absent, yet it orders. On the 
other hand, Tb2Ti207 and Ho2Ti207 in which crystal 
fields and anisotropy are clearly important, do not order 
and it is likely that these perturbations in fact inhibit 
the occurrence of long range order by competing with 
important interactions other ^M than nearest-neighbor 
Heisenberg antiferromagnetic exchange ||2y,^,^. Some 
other interesting recent results have been found in the 
Yb2Ti207 and Dy2Ti2 07 pyrochlores [||. 

It is also useful to compare Gd2Ti2 07 with other 
Heisenberg pyrochlores such as Y2M02O7 ]2^-|24[] and 
Heisenberg spinels such as ZnFe204 [£l[. Y2M02O7 is 
a well-known geometrically frustrated antiferromagnet 
spin glass material with Ocw /Tf > 10, Tf being the spin 
freezing temperature of 21K pJ. Here too, only a spec- 
ulation can be offered for the differences as follows. Be- 
cause of the high level of degeneracy across the zone, the 
nearest-neighbor Heisenberg pyrochlore antiferromagnet 



is expected to be fragile against a small, random disorder 
level, X, and will have a propensity to develop a disorder- 
driven spin glass ground state the smaller the perturba- 
tions {H'} beyond the nearest-neighbor exchange inter- 
action is |n,pO| . We expect that the critical disorder level 
for the Neel to spin-glass transition, Xc, will go to zero 
as {H'} goes to zero m. For example, in Y2M02O7, 
there is preliminary evidence that the second neighbor 
exchange parameter, J2, is only a few percent of Ji [ p2[ . 
In Gd2Ti207 on the other hand the leading corrections 
{H'} are dipolar interactions. Dm, and of order 20% of 
Ji. In other words {H'}/Ji is not smaU in Gd2Ti207 
and the anisotropy of the dipolar interactions will possi- 
bly introduce sizeable stabilizing anisotropy gaps to the 
spin-wave excitations out of the selected long range or- 
dered ground state. Both the relative size of Ddd/Ji 
and the "spin-holding" effect of the anisotropy of dipo- 
lar interactions will result in a much increased Xc com- 
pared to more isotropic Heisenberg systems with small 
{H'} . In summary in this picture weak disorder drives 
the spin glass transition in Y2M02O7 while the strong 
and anisotropic dipolar interactions "helps" stabilize long 
range order in Gd2Ti207. In this context the existence 
of a very weakly dispersive line along [111] restored by 
order-by-disorder, or perturbative J2 and J3, would sug- 
gest that, as in the fee antiferromagnet [|l6[, weak random 
disorder would rapidly drive Gd2Ti207 into a spin glass 
state. 

Finally, it is interesting to compare the behavior of 
Gd2Ti207 with the frustrated ZnFe204 antiferromagnet 
spinel where Fe'^^ is a ^6*5/2 closed shell ion for which 
single-ion anisotropy should be negligible as is the case 
for the *S'7/2 Gd^+ ion in Gd2Ti207 . In the insulating 
normal Heisenberg spinel ZnFe204, where the Fe^+ mag- 
netic moments occupy a lattice of corner-sharing tetra- 
hedra, muon spin relaxation and neutron studies have 
revealed that long range antiferromagnetic order (LRO) 
develops below Tn — 10. 5K. However, already at tem- 
peratures of about T « IOTn a short range antiferro- 
magnetic order (SRO) develops which extends through 
w70% of the sample volume just above Tn- Below Tn 
antiferromagnetic SRO and LRO coexist. At 4.2K still 
w 20% of the sample are short range ordered. The re- 
gions exhibiting SRO are very small « 30A. The phys- 
ical origin of the SRO as well as partial glassy behavior 
in ZnFe204 remains an enigma. Hence, while it appears 
that Gd2Ti207 displays conventional antiferromagnetic 
long range order and Y2M02O7 shows full-blown spin- 
glass behavior, ZnFe2 04 exhibits a combination of both 
short and long range antiferromagnetic order, as well as 
spin-glassy behavior. The origin of the difference between 
the Gd and Fe based pyrochlore lattice antiferromagnets 
in terms of their coexistence of long range order and spin- 
glassy behavior is not known. Possibly different range 
of interactions, presence of strong dipolar anisotropy in 
Gd2Ti2 07 compared to a much more overall isotropic 
spin-spin interaction in ZnFe204 may play some role. In 



light of this, it would be interesting to study in further 
detail the magnetic properties of Gd2Ti207 using muon 
spin relaxation and neutron scattering methods. 



VI. CONCLUSION 

Evidence has been presented from a.c. and d.c. suscp- 
tibility and specific heat measurements, that the frus- 
trated Gd2Ti207 insulating pyrochlore exhibits a transi- 
tion to a long range ordered state at 0.97K as opposed to 
a spin glass or spin liquid state as often observed in other 
pyrochlore materials. From specific heat measurements, 
short range magnetic correlations have been found to ex- 
tend to T > 30Tc and the entropy removal below Tc is 
only about 50%. From a mean field theoretical study it 
is concluded that no long range order should exist for the 
pyrochlore lattice for nearest neighbor antifcrromagnetic 
interactions, Ji, only, even upon inclusion of long range, 
anisotropic dipolar couplings, Ddd- Long range order at 
various commensurate or incommensurate wave vectors 
is predicted to occur only upon including a finite second, 
J2, and/or third, J3, nearest- neighbor exchange interac- 
tions beyond Ji and Ddd- Long range order could also 
be driven by thermal and/or quantum fluctuations via 
an order-by-disorder mechanism. The exact wave vector 
depends on the relative signs and magnitudes of J2, J3 
and Ddd- It would be of interest to investigate further 
the nature of the ordered state in zero and applied fields 
in Gd2Ti207 by neutron scattering and muon spin relax- 
ation methods. Finally, we argued above that the related 
fee antiferromagnet material, cubic-Gd203, should be re- 
considered as a geometrically frustrated antiferromagnet 
and is worthy of further study. 
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TABLE I 

Presence or Absence of Long Range Order in RE2Ti2 07 
Pyrochlores. 
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FIG. 1. Pyrochlore lattice of corner-sharing tetrahedra. 
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FIG. 2. (a) Inverse molar susceptibility, 1/x, of Gd2Ti2 07 
against temperature in the temperature range T — [2 — 25] 
K, and in the temperature range T = [2 — 300] K in the inset, 
(b) Inverse molar susceptibility, l/x-, of (Gdo.o2Yo.98)2Ti207 
against temperature in the temperature range T = [2 — 25] K, 
and in the temperature range T = [2 — 300] K in the inset. 
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FIG. 5. Magnetic specific heat, Cm, (obtained by subtract- 
ing Ci from Cp) against temperature. The solid line repre- 
sents the theoretical fit (see text for details). The inset shows 
a blown-up region of Cm in the low-temperature regime. 
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FIG. 3. Real part of the ac susceptibility, x', vs. tempera- 
ture measured at diflerent frequencies. 
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FIG. 4. Specific heat, Cp, of Y2Ti2 07 as function of tem- 
perature. The solid line corresponds to the lattice spe- 
cific heat, Ci, estimated from the measurements on the 
non-magnetic Y2Ti207. 
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FIG. 6. Largest eigenvalue Amax(q) as a function of 
wavevector q for Ddd/ Ji ~ 0.2 and J2 — Js = 0. A de- 
generation line occurs for q in the star of the [111] direction. 
The small "ripples" seen on the degeneration lines along (111) 
and (111) directions are due to the finite number (500) near- 
est-neighbors considered in the dipolar interactions. When 
considering more than 10 nearest-neighbors, the maximum of 
A(q) always occurs on the star of [111] with the amplitude of 
the modulations due to the dipolar cut-off continuously de- 
creasing as the number of nearest-neighbors is increased to 
infinity. 
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